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ABSTRACT 

We present Chandra and HST observations of the ultraluminous X-ray source (ULX) IC 342 X-l. 
The Chandra and HST images are aligned using two X-ray emitting foreground stars. The astrometry 
corrected position for X-l is RA. = 03 h 45 m 55!61, Decl. = +68°04'55'.'3 (J2000) with an error circle of 
0.2". One extended optical source is found in the error circle, which could be the optical counterpart 
of X-l. The source shows an extended feature in HST images at long wavelengths, which is likely 
to be a superposition of two point sources, although it is possible that the dimmer one could be a 
jet. Both sources are much redder than typical for ULX optical counterparts. The brighter one has 
an absolute magnitude My — —5.2 ± 0.2 and (B — V)q = 0.66 ± 0.13 and the dimmer star is not 
detected in B and has (B — V)o > 2.1. Their colors are consistent with an F8 to GO lb supergiant 
or a carbon star, respectively. However, it is likely that part or most of the optical emission may be 
due to X-rays reprocessed by the companion star or the accretion disk. The stellar neighborhood of 
IC 342 X-l lacks O stars and has a minimum age of ~ 10 Myr. This excludes the possibility that the 
surrounding nebula is powered by an energetic explosion of a single massive star that formed a black 
hole. We suggest that the nebula is most likely powered by an outflow from the X-ray source. 
Subject headings: black hole physics — accretion, accretion disks — galaxies: individual (IC 342) — 
galaxies: stellar content 



1. INTRODUCTION 

Ultraluminous X-ray sources, the brightest nonnu- 
clear X-ray point sources in external galaxies, still hide 
their nature from us. Their luminosity of 3 x 10 39 
to 10 41 ergs s _1 , inferred assuming isotropic emis- 
sion, is too high to be powered by accretion onto 
stellar mass black holes that are formed in the col- 
lapse of normal stars with normal metallicity, and thus 
may suggest the presence of intermediate mass black 
holes (IMBHs: IColbert fc MushotzkvlH999 : Kaare t et all 
2001 ). However, if the em ission is beamed ( King et al.l 
20011: lK5rding et all |2002|) or truly super-Eddington 
(IWatarai et all 120011: fBegelmarJkOO^ . IMBHs are not 



required. 

Multiwavclcngth observations have provided pow- 
erful means to dete rmine the nature of ULXs. 
iPakull fc Mirionll (j2003j) carried out an optical survey of 
bright ULXs and discovered that several ULXs are as- 
sociated with optical nebulae. These nebulae have very 
large spatial scales of a few hundred parsecs, and present 
line emission. In particular, He n A4686 emission line 
was detected in the nebula surrounding Holmberg II X-l, 
and was used to determine that the flux of reprocessed X- 
rays from the centra l ULX was at least 4 x 10 39 ergs s _1 
(|Kaaret et alJl2004h . Identification of an optical coun- 
terpart of a ULX has been accomplished in several cases. 
However, it is hard to determine whether or not the op- 
tical light from the companion star is contaminated by 
the outer region of the accretion disk, thus stellar classi- 
fication based on optical colors may not be reliable. The 
far- ultraviolet spectrum of NGC 5204 X-l suggests a B0 
lb companion, b ut contains some peculiar emission lines 
(|Liu et al.l 2004). Identification of a giant branch com- 
panion has been made for the ULX X41.4+60 in M82 
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based on the detection of a 62 day X-ray modulati on 
interpreted as the or bital period ([Kaaret et al.ll2006aH bl; 
iKaaret fc Fengl l2007f ) . Identification of the companion 
stars for more ULXs would enable us to place constraints 
on the binary evolution of these systems. 

IC 432 is a nearly face-on starburst spiral 
galaxy. IC 342 X-l was discovered with Einstein 
([Fabbinao fc Trinch ieri 19815) and confirmed as a poin t 
source with ROSAT observations (iBreeman et al.lfl99 
Using ASCA observations in 1993, Okada et alj (|199< 
found significant variability on time scales of hours from 
X-l, suggesting it was a compact object. The source 
changed its spectrum dramati cally in a 2000 ASCA 
observation (|Kubota et al.|[200Th . similar to the spectral 
state transitions found in Galactic black hole binaries. 
A possible X-ray periodicity of 31 or 41 hours from 
X-l was also detected in t he same ASCA observation 
in 2000 (jSugiho et al.ll2001l ). X-l is associated with a 
"tooth" -shaped optical nebula (jPakull fc Mirionil [2003] ) 
with e mission lines simi l ar to a supernova remnant 
(SNR: IPakull & : Mirionl [200l IRoberts et all [200l 
lAbolmasov et al.ll2007t ). iRoberts et al.l (|2003h estimated 
that the initial energy required to power the nebula was 
at le ast 2 x 10 51 ergs, gr eater than a canonical SNR. 
Also, IRoberts et al.l (|2003l ) found two regions of [O in] 
emission located oppositely on the shell of the nebula, 
which could be X-ray photoionized and t hus indicate 
anisot ropic emission from X-l. Instead, iGrise et al.l 
(|2006| ) found that the [O III] emission closely follows the 
Ha emission and is smoothly distributed in the n ebula , 
implying it is shock excited. lAbolmasov et all (|2007f ) 
found evidence for photoionization in the central regions 
of the nebula. 

To search for an optical counterpart of IC 342 X-l, 
we obtained Chandra and HST observations. The ob- 
servations and data are described in § [2j The Chandra 
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and HS T images are aligned by matching two foreground 
stars in both images (§ [3]) , which leads to the identifica- 
tion of a unique optical counterpart of X-l (§ SJ. The 
photometry for the optical counterpart of X-l and its 
nearby starforming region is performed in § A multi- 
wavelength spectrum from optical to X-ray is presented 
in §[B] An emission line image obtained using HST of the 
tooth nebula is shown in § [71 All results are discussed 
in § M We adopt a distance to IC 342 of 3.3 ± 0.3 Mpc, 
which is es timated from obs ervations of 20 Cepheids in 
the galaxy (|Saha et al.ll2002D . 

2. OBSERVATIONS AND ANALYSIS 

HST observed the IC 342 X-l region with the Ad- 
vanced Camera for Surveys (ACS) under the proposal 
10768 (PI: P. Kaaret) in December 2005. Four observa- 
tions, with an exposure of 920 s each, were made using 
the Wide Field Camera (WFC) with a broadband filter 
F625W (Sloan Digital Sky Survey R). Each observation 
had the same orientation, but had IC 342 X-l located in 
a different corner of the image in order to create a large 
mosaic image around X-l. Another four WFC images 
were obtained with X-l at the center using three broad- 
band filters, F435W (Johnson B), F555W (Johnson V), 
F814W (Broad I), and a narrowband filter, F658N (Ha + 
[N n] ) . Each of the three broadband images has an expo- 
sure of 1800 s, and the narrowband one has an exposure 
of 2400 s. The fields of view for all images are shown 
on an STScI Digitized Sky Survey (DSS) blue image of 
IC 342 (Figure [T] top). We used the standard drizzled 
ACS images for all analysis. 

The Chandra X-ray Observatory observed IC 342 X- 
1 on 2006 Nov 12 (ObsID 7069; PI: P. Kaaret) with 
the Advanced CCD Imaging Spectrometer (ACIS). An 
X-ray image was created in the 0.3-8 keV range, cover- 
ing the field of view of the mosaic F625W image. We 
performed source detection using the wavdetect tool in 
CI AO 3.3.0.1 applying an exposure map calculated for 
a power- law spectrum with a photon index r = 1.5 and 
Galactic absorption iVn = 3.02 x 10 21 cm~ 2 . 

3. ASTROMETRY CORRECTION 

Each HST im age was aligned to t he sources in the 
2MASS catalog (|Skrutskie et al.ll2006l ) using the Graph- 
ical Astronomy and Image Analysis Tool (GAIA) to cor- 
rect the absolute astrometry. We selected catalog stars 
that appear as pointlike sources in an HST image, and 
found 22 to 35 counterparts in each image. The rms 
(root mean square) deviation of the HST source posi- 
tions versus the 2MASS catalog positions is <2 pixels 
(w 0.1"). The astrometric accuracy for 2MASS stars 
with 13 < Ks < 15 (ones used here) has an rms of 
around 0.1" along t he right ascension (see Fig. 20 in 
ISkrutskie et al.ll2~006h . which corresponds to a 90% posi- 
tional error of 0.22" (assuming a Rayleigh distribution). 
After correction, the statistical accuracy of the alignment 
of each HST image relative to the 2MASS frame is about 
<5 H = 0.22"/V^V w 0.04" at 90% confidence, where N (= 
22 to 35) is the number of counterparts. 

The four F625W images, after astrometry correction, 
are combined into a mosaic image (see Figure [TJ bottom) 
using the IRAF package MSCRED. A flux offset was im- 
posed on each image in order to equalize the mean flux in 
the overlapping image regions. This helps remove image 



artifacts at the image boundaries. 

Six counterparts are found between the X-ray image 
and the F625W mosaic image. Two of them are located 
on a different ACIS CCD chip from X-l, thus not suit- 
able for astrometry. Another two sources have only 15 
and 17 photons detected by Chandra, and the bright- 
est pixel around the source region contains only 2 and 
4 photons, respectively; we decided not to use them for 
astrometry due to their low significance. The other two 
sources have 617 and 74 photons detected by Chandra, 
respectively, with a significance of 177cr and 18ct given 
by wavdetect. These two sources are used for astrome- 
try, and their HST and Chandra positions are shown in 
Table [TJ Source 1 is highly saturated in the HST image, 
and source 2 is moderately saturated. Their HST po- 
sitions are obtained by measuring the crosspoint of the 
diffraction spikes, and the errors are estimated to be less 
than 0.05". The X-ray positions and statistical errors are 
obtained from wavdetect. There are systematic offsets 
in the positions reported by wavdetect caused by the 
shape of the point spread function. To estimate these, 
we performed simulations using MARX of sources lo- 
cated at the observed positions, fit the simulated images 
using wavdetect, and calculated the offset between the 
sky position input to the simulation and the position re- 
ported by wavdetect. These systematic offsets are listed 
in column (7) and (8) of Table [TJ for each source. The 
corrected Chandra positions can be obtained by adding 
the positions quoted in Table [TJ and the c orresponding 
offset. Source 1 is referred as source 1 1 in I Bauer et al.l 
(|2003h and X12 in lKond (|2003h . It appears to be a fore- 
ground star. Source 2 appears in several catalogs and 
it s colors and X-ray to o ptical ratio (> 0.001; as defined 
in lMaccacaro et al"1ll982[ ) are also consistent with a fore- 
ground star. 

From wavdetect, the position of X-l is obtained as 
R.A. = 03 h 45 m 55?636, Deck = +68°04'55"41 (J2000) 
with a statistical error of <5xi = 0.02" (90% confidence). 
The system offset in the wavdetect position estimated 
using a MARX simulation as described above is Ara = 
0.03", Adcc = 0.01". The absolute astrometric error for 
Chandra is 5c = 0.6" at 90% confidence 2 . 

We used the geomap task in the IRAF package IM- 
MATCH to align the X-ray image to the mosaic F625W 
image by fitting the positions of the two sources in Ta- 
ble [TJ and then applying the best-fit transform parame- 
ters onto X-l using the task geoxytran. This gives a best 
position of R.A. = 03 h 45 m 55 i !612, Deck = +68°04'55 / .'29 
(J2000) for X-l. The uncertainty in the fit using the 
two sources is Sat = \J &\ + <5 2 /2 = 0.21", where 5i and 
82 are the positional errors for sources 1 and 2 in Ta- 
ble [TJ Taking into account the error of X-l, we esti- 
mate the final 90% error radius on the position to be 
\J <5fr t + <$xi = 0.21". We note that a consistent position 
can be obtained for X-l by matching source 1 only, but 
this yields comparable positional error because of the un- 
certainty in the rotation angle between the images. The 
F625W image at a region around X-l is presented in Fig- 
ure [21 The uncorrected Chandra error circle is displayed 
as a dashed circle, and the corrected position of X-l with 
the error circle is shown as a small solid circle. 

2 http : //exc . harvard . edu/cal/ASPECT/celmon/ 
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Fig. 1. — Top: DSS blue band image of IC 342. The large rectangle indicates the field of view of the mosaic F625W image, and the 
small one indicates the field of view for the F435W, F555W, F814W, and F658N images. The arrow points north and has a length of 2'. 
Bottom: Mosaic F625W image. The positions of X-l and the two sources for astrometry are labeled. The arrow points north and has a 
length of T. 



TABLE 1 

HST AND Chandra POSITIONS OF THE TWO SOURCES FOR. ASTROMETRY CORRECTION. 





HST 




Chandra 








No. 


R.A. Decl. 


R.A. 


Decl. 


5 


A RA 






(J2000) (J2000) 


(J2000) 


(J2000) 


(") 


(") 


(") 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


03 46 06.592 +68 07 05.42 


03 46 06.627 


+68 07 05.49 


0.09 


0.00 


0.07 


2 


03 46 25.985 +68 04 19.87 


03 46 26.004 


+68 04 19.94 


0.41 


0.24 


0.06 



Note. — <5 is the statistical 90% error radius for Chandra positions. A is the systematic 
offset of the wavdetect position caused by the asymmetry of the point spread function. 
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Fig. 2. — Images in the region around IC 342 X-l. In the F625W image, the dashed circle indicates the directly measured Chandra 
position of X-l with an uncertainty of 0.6", the solid circle indicates the corrected X-l position with an uncertainty of 0.21", and the two 
possible optical counterparts, stars A and B, are labeled. The arrows points north and are 1" long. 



F435W 





Pixel Pixel 

Fig. 3. — Profiles of brightness along the X (solid line) and Y (dotted line) directions of the CCD across the brightest pixel of star A. X 
and Y profiles are centered at their best-fit centroids. There is a secondary component on the X profile at long wavelengths, corresponding 
to star B in Figure [2] 
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TABLE 2 

Magnitudes and colors of the two stars near IC 342 X-l 





Star A 


Star B 


Filter/Band 


m mo 
(mag) (mag) 


m mo 
(mag) (mag) 



ST magnitudes 



F435W 24.68 ± 0.07 22.40 ± 0.09 > 28.1 > 25.8 

F555W 24.15 ± 0.06 22.33 ± 0.07 26.08 ± 0.17 24.30 ± 0.18 

F625W 23.95 ± 0.06 22.48 ± 0.07 25.44 ± 0.13 24.02 ± 0.15 

F814W 23.77 ± 0.05 22.68 ± 0.05 24.77 ± 0.08 23.73 ± 0.09 



VEGA magnitudes 



B 


25.29 ± 0.10 


23.01 ± 0.10 


> 28.6 


> 26.3 


V 


24.17 ± 0.08 


22.35 ± 0.08 


26.0 ± 0.2 


24.2 ± 0.2 


R 


23.45 ± 0.07 


21.98 ± 0.07 


24.6 ± 0.4 


23.2 ± 0.4 


I 


22.50 ± 0.06 


21.41 ± 0.06 


23.62 ± 0.18 


22.58 ± 0.18 


M v 




-5.2 ± 0.2 




-3.4 ± 0.3 


(B-V) 




0.66 ± 0.13 




> 2.1 


(V-R)o 




0.37 ± 0.11 




1.0 ± 0.4 


(R-I)o 




0.56 ± 0.10 




0.6 ± 0.4 



Note. — m is the observed magnitude and mo is the Galactic reddening 
corrected magnitude assuming E(B — V) = 0.558. Magnitudes are calculated 
based on the spectrum of an F8 I supergiant for star A, and of an M5 III gaint 
for star B. 



4. OPTICAL COUNTERPART OF X-l 

After the astrometry correction, a likely optical coun- 
terpart falls into the error circle of X-l, see Figure El 
The source is larger than the point spread function of 
the ACS with the F625W filter. Figure [3] shows the pro- 
files of brightness along the X and Y pixel directions 
passing through the brightest pixel in the source for all 
broadband images. It is obvious that there is an excess 
along the X profile which is asymmetric and more ex- 
tended than the Y profile in all the images except the 
F435W image. This excess indicates that the source is 
either truly extended, a star with an associated extended 
source, or two stars. We consider the last possibility most 
likely and label the bright core of the source lying well 
within the error circle of X-l as "star A" and the other 
feature as "star B" in Figure [H Star B is hardly seen 
in the F435W image, but becomes more and more sig- 
nificant at longer wavelength. It lies about 0.16" away 
from star A on the sky plane, just on the edge of the 90% 
confidence error circle of X-l. 

5. PHOTOMETRY 

Simple aperture photometry is not adequate to mea- 
sure the magnitudes of stars A and B since they are not 
fully resolved. Instead, we performed a simultaneous fit 
to the two sources with a point spread function, measured 
from bright isolated stars in the same image, to obtain 
their flux. The fit was performed in the Interactive Data 
Language (IDL) using an algorithm adapted from the 
IRAF package DAOPHOT by W. B. Landsman. Star B 
is required for an adequate fit except for the F435W im- 
age. The flux upper limit of star B in the F435W image 
was estimated from the total residuals around its position 
after star A was subtracted. Magnitudes are calculated 
from the flux measurements using the SYNPHOT pack- 
age, assuming an F8 I supergiant spectrum for star A 
and an M5 III giant spectrum for star B, respectively. 
The Galactic extinction along the line of sight to IC 342 
is E(B — V) = 0.558, which i s calculated from du st maps 
derived from COBE data (jSchlegel et al.lll998[ ). The 



Galactic extinction used here is consistent with the total 
ex tinction to IC 342 of E(B-V) = 0.63± 0.10 estimated 
bv lSaha et all (12002ft . iGrise et alj (|2006( ) estimated a lo- 
cal extinction of E(B — V) = 0.26 from the H a /Hp ratio 
in the n ebula, indicating t hat m ost extinction is from our 
Galaxy. iFingerhut et al.l (|2007ft also found that the ma- 
jority of the extinction along the line of sight to IC 342 
was Galactic. Therefore, we choose to not consider extra- 
galactic extinction in calculating the reddening corrected 
colors. 

Observed and reddening corrected ST magnitudes and 
VEGA magnitudes are listed in Table [H To estimate 
systematic errors caused by different spectral types and 
band translation, we calculated magnitudes using spectra 
of different stars and added the errors linearly onto the 
results as systematic errors. Spectral types used to esti- 
mate the systematic errors include F2 I, GO I, and G2 I 
for star A, and all M III for star B. The absolute V mag- 
nitudes are calcul ated assuming (m — M) = 27.58±0.18 
dSaha et all l2"00l) . Colors of (B - V)q, (V - R) a , and 
(R — I)q for both stars are also listed in Table [21 Both 
star A and star B appear to have constant flux in the 
four F625W images, which were obtained over an inter- 
val of nearly 5 hours. The V magnitude and colors of 
star A are consistent with an F8 to GO lb supergiant, 
except for some excess brightness in the / band. The 
magnitude and colors of star B indicates it could be an 
M II giant; its red color suggests that it is similar to a 
carbon star ()Mendoza k. Johnsonlll965l ). 

A true color image of the region within 20" of X- 
1 is presented in Figure |4] composed of the F814W, 
F555W and F435W images, respectively, for the red, 
green and blue channels. Using the I DL DAOPHOT 
and following the procedures described in Kaaret ( 2005) , 
we present a color-magnitude diagram of Mpsssw ver- 
sus Mp435\v — A/F555W in Figure \E\ for bright stars in 
the region that shows in Figure [U To exclude extended 
sources, image defects and foreground stars, we selected 
sources that had peak values less than 0.75 and a local 
fit \ 2 < 3. The magnitudes are extinction corrected but 
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Fig. 4. — True color image of a 40" X 40" region around IC 342 X-l. The black belt near the bottom is caused by the CCD gap. The 
F814W, F555W, and F435W images are used for the red, green, and blue channels, respectively. The arrow points north and has a length 
of 5". 

there might be additional reddening in the host galaxy 
or internal to the source, which is not considered. Typ- 
ical errors are 0.2 for magnitudes and 0.16 for colors. 
Stars with absolute magnitudes brighter than —8 and 
also stars with colors bluer than —0.4 are likely fore- 
ground stars. Isochrones for stellar populations with a 
metallicity of Z = 0.019 at various ages (Girardi 2006, 
|http : //pleiadi . pd . astro . it ) are also shown in Fig- 
ure [5] 

6. MULTIWAVELENGTH SPECTRUM 

In the Chandra observation, there are 0.465 counts per 
frame detected on a 3 x 3 pixel island around X- 1 . Follow- 
ing The Chandra ABC Guide to Pileup 3 , this detected 
rate requires an input rate of 0.92 or 2.00 counts per 
frame assuming a grade migration probability a = 0.5, 
corresponding to a count lost of 49% or 77%, respec- 
tively. Using PIMMS, we find that the observed fluxes 
from ASCA and XMM are roughly consistent with the 
detected count rate. However, the spectral shape and 
the luminosity state cannot be determined. 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

M F435w - M F55sw [VEGAmag] 

Fig. 5. — Color-magnitude diagram for stars shown in Figure [4] 
Star A is indicated with a star symbol. Isochrones of 10, 20, and 
100 Myr stars are superimposed. We note that the magnitudes are 
found in the HST filter bands but quoted in the Vega magnitude 
system. 



li http : //asc . harvard . edu/clao/download/doc/pileup_abc . ps 
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Fig. 6. — Multiwavelength spectrum of IC 342 X-l. Corsses 
and squares are reddening corrected optical fluxes for star A and B 
respectively, and lines are absorption corrected best-fit models for 
the X-ray spec tra. The solid l ines r epresent the two ASCA spectra 
obtained from Kubota ct al. (2001) and the the dashed line is the 
XMM spectrum from lBauer et al.l 12003) . We note that the ASCA 
spectra could be contamin ated by surrounding sources by a factor 
of 20% HBauer et alj|2003h . 

IC 3 42 X-l is a highly variable source. iKubota et al.l 
(2001) reported that the source changed its X-ray spec- 
tral shape dramatically from 1993 to 2000 depending 
upon two ASCA observations. XMM observed the source 
in 2001 and found the source spectrum was similar to 
that observed in the 2000 ASCA observation. A mul- 
tiwavelength spectrum consisting of optical and X-ray 
data is presented in Figure O The optical spectra show 
a form of S u oc ^-^±0.12 for star A and ^ K ^- 3 .2±o.4 

for star B. We note that the ASCA spectra could be con- 
taminated by surrounding sources and the inferred flux 
could be overestimated by a factor of 20% (| Bauer et al.l 
12001 . 

The X-ray to optical flux ratio, defined as 
log(/x//y) = log/x + rny/2.5 + 5.37, where fx is 
the 0.3-3.5 keV flu x in ergs" 1 cm" 2 s" 1 a nd my is the 
visual magnitude (Maccac aro et al.lll982D . is found to 
be 3.0-3.7 assuming star A is the counterpart or 3.8-4.5 
for star B, given different X-ray spectra and the HST 
my. This is larger than the ratios for normal stars, 
normal galaxies, clusters of galaxies, AGNs, and BL 
Lac objects found in the Einstein Extended Medium 
Sensi tivity Survey (jMaccacaro et alj Il988t IStocke et al.1 
1991). Similarly, the two-point optical to X-ray spectral 
index computed between the V band and 1 keV is 
aox = 0.2 - 0.4 for star A and a OX = -0.05 - 0.14 for 
star B, which are flatter than the values foun d in those 
objects ()Stocke et alj|1991t Eajidt et al.ll2001l) . 

To compare with Galactic X-ray binaries, we com- 
puted the X-ray to optical flux ratio defined in the form 
£ = B + 2.5 log Fx, where B is the reddening cor- 
rected B magnitude and Fx is the 2-10 keV observed 
X-ray flux in /iJy. We obtain £ = 20.4 — 22.1 for star 
A and X-l using its X-ray spectra at different emission 
states. Since star B is not detected in the B band with 
these HST observations, we estimate the £ for star B is 
greater than 24. These numbers are greater than that 
from most high-mass X-ray binaries, but well consistent 
with the range of £ found in low -mass X-ray binaries 
van Paradiis fc McClintocldH995l) . 



The multiwavelength spectrum of X-l, no matter 
whether star A or star B is assumed to be its optical 
counterpart, is unlike the spectra of quasars or blazars. 
The spectrum is consistent with an X-ray binary in IC 
342 X-l. 

7. LINE EMISSION FROM THE NEBULA AND X-l 

Figure [7] shows the F658N image in the region around 
IC 342 X-l. This filter includes lines from [N 11] at A6583 
and A6548 as well as Ha at A6563, and their relative 
transparency is 1.0, 0.5 and 1.0, respectively. The error 
circle of X-l is shown and the optical counterpart of X-l, 
discussed above, is clearly visible in the image. We note 
that the optical counterpart of IC 342 X-l lies very close 
to the geometric center of the emission line nebula. This 
strengthens the association between the nebula and the 
X-ray source and suggests that the X-ray source has not 
moved significantly since the creation of the nebula. 

The nebula has an angular scale of around 6" in diam- 
eter, corresponding to 100 pc at a distance of 3.3 Mpc. 
We estimated the total flux of the nebula using an ir- 
regular shaped region encompassing the regions bright 
in the F658N image but with stars apparent in the con- 
tinuum image (F625W) removed. We used the F625W 
image to estimate the continuum. Since the line and con- 
tinuum are shown in both filters, we modeled the emis- 
sion as a flat continuum plus line emission and solved 
algebraically for the continuum and line strength given 
the count rates in the two filters. We estimate that 
the reddening corrected line flux from the nebula is 
(4.7 ±0.3) x 10~ 14 ergs" 1 cm -2 s _1 . This flux is consis- 
tent with the sum of th e Ha and [N 11] fluxes found by 
lAbolmasov et all ((2007). We note that this flux i s signif - 
icantly higher than that quoted bv lRoberts et alj (2003), 
but we find that our flux before the red dening correc- 
tion i s consistent with that quoted by iRoberts et al.l 
(120031) as redden i ng co rrect ed. Using the line r atios in 
lAbolmasov et al.l(|2007t ) and lRoberts et all (|2003l ). we es- 
timate that the Ha A6563 flux from the nebula is about 
(2.5 ± 0.3) x 10~ 14 ergs" 1 cm -2 s -1 . Using the same 
continuum subtraction procedure, we do not detect any 
Ha line emission from the optical counterpart of X-l and 
place an upper limit on the reddening corrected flux of 
1.0 x 10 -16 ergs -1 cm -2 s -1 . 

8. DISCUSSION 

The stellar population within 300 pc (20") of X-l ap- 
pears to contain no O stars. The color-magnitude dia- 
gram (Figure [5]), suggests that the minimum stellar age 
is ~10 Myr if the metallicity is close to Z = 0.019. Com- 
paring to isochrones at a lower metallicity of Z = 0.004, 
the minimum stellar age is about 20 Myr. The region ap- 
pears to contain relatively young stars with ages of tens 
of Myr and also an older population with ages of several 
Gyr. 

IC 342 X-l appears at the center of the nebula (Fig- 
ure [7]), which strongly suggests an association. The rel- 
atively old age of the stellar population in the vicinity 
of the nebula makes it unlikely that the nebula is pow- 
ered ( even in part) by ioniza tion from one or more O 
stars. lAbolmasov et al.1 (|2007|) suggested that the shock 
velocity in the nebula is 20 to 100 km/s, implying the 
nebula is younger than 2 Myr. This age is an overesti- 
mate since the expansion would be slowing down. This 
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Fig. 7.— F658N image of the region around IC 342 X-l. This filter includes lines Ha A6563, [N II] A6583, and [N II] A6548. The error 
circle for X-l is shown. The arrow points north and has a length of 1". 



suggests that if the nebula was created in a supernova 
event, the progenitor star must has a lifetime close to 
10 Myr. The high luminosity of X-l suggests it contains 
a black hole. To form a black hole in an energetic explo- 
sion, the progenitor must b e very massive (e.g., at least 
40 M & in lHeger et al.|[2003h and thus has a lifetime sig- 
nificantly shorter than 10 Myr. Therefore, it is unlikely 
that the nebula was powered by an unusually powerful 
explosion in which a black hole was formed. 

A more plausible alternative is that the nebula is pow- 
ered by an outflow from the accreting compact object i n 
the ULX (|Pakull fc Mirionil l200l iRoberts et aflf2003l) . 
The nebula appears to be predominantly shock powered 
with a shock velocity in the range 20-100 km/s and a 
total power of 3 x 10 39 ergs s , but with high excita- 
tion lin es possibly powered by photoionization near the 
center (jAbolmasov et all 12007( 1. This power is 20% to 
50% of X-ray lumi nosity, depending on the X-ray state 
(Oka da~et alJ ll998). and is reasonable in that X-ray bi- 
naries appear to have power in their outflows which is 
comp arable to or may even exceed their X-ray emis- 
sion (jGallo et alJl2005h . The emission line morphology 
in Figure [7] indicates that the nebula is not a complete 
shell. The morpholo gy is somewhat simil ar to that of the 
W50/SS433 system |Boumis et all 12007( 1 and may sug- 
gest a jet orientation along an axis slightly west of north. 
Given the strong radio emission from W50/SS433, radio 
observations of IC 342 X-l would be of great interest. 

We find one optical source in the error circle of IC 342 
X-l. The source appears extended at long wavelengths 
and is either two stars or a star with an associated feature 
such as a jet. We consider the relativistic jet interpreta- 



tion unlikely because their spectra are quite steep, unlike 
the o ptical relativistic jets of black hole such as those of 
M87 (|Perlman et al.ll2001f ). The projected separation in 
IC 342 would be 2.6 pc, and the true separation may be 
larger. This is too large to be an optical compact jet 
found in the hard state of stellar mass black ho les, which 
usually has a spatial scale of AU (jFenderll2006( ). We find 
interpretation in terms of two superimposed stars more 
likely. The brighter star (star A) lies well within the 
error circle for IC 342 X-l. The dimmer star (star B) 
lies at the edge of the 90% confidence error circle of X-l 
and could, instead, be the counterpart of X-l. Allowing 
either star as the counterpart, the X-ray to optical flux 
ratio indicates that the source is not a background blazar 
or quasar, but is consistent with an X-ray binary in IC 
342. 

If the companion star dominates the optical emission 
despite the high X-ray flux, then the spectral types of 
the stars can be determined. The colors and absolute 
magnitude of star A are consistent with classification as 
an F8 to GO lb supergiant, which has a mass of about 
10 Mq. There is some excess I band emission, but this is 
often seen in supergiants and is a signature of a circum- 
stellar disk or a stellar wind. With this classification, the 
age of star A would be roughly 50 Myr. The evolution- 
ary status of star A, on the supergiant branch, would be 
compatible with the large mass transfer rates needed to 
power the inferred luminosity of the ULX. 

The colors and absolute magnitude of star B are con- 
sistent with classification as a carbon star. The progen- 
itor mass of such a star is likely of intermediate mass, 
1 — 5Mq, and, thus, would require a rather old binary 
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system. The carbon star phase is short-lived, < 1 Myr. 
High lumin osity carbon stars have high mass loss rates 
(|Juralll99ll ) which would be sufficient to power the ULX. 

However, the optical emission could also be dominated 
by reprocessed X-rays on the companion star or on the 
accretion disk, or by the intrinsic emission of the disk. In 
these cases, spectral classification of the companion star 
via its colors is not reliable. We discuss these possibilities 
in following. 

Assuming the X-ray emission is isotropic, the X-ray 
luminosity intercepted by the companion is 



47ra 2 



Lx , 



(1) 



where R2 is the radius of the companion star, a is the 
binary separation, and Lx is the total X-ray luminos- 
ity. The high luminosity of ULXs indicate that they ac- 
cr ete via Roche-lob e overflow. Substituting equation (2) 
in lEggletonl (JT983) in which R^ja is determined by the 
companion to compact object mass ratio q to an accuracy 
better than 1% for any q, we have 



L r = 



Lx 
4 



0.49 g 2 / 3 
0.6(? 2 / 3 + ln(l -f 



„1/3A 



(2) 



A plot of L c versus q at different X-ray luminosities is 
presented in Figure [5J We can see that for a source like 
IC 342 X-l, which has Lx ~ 10 40 ergs s _1 , the compan- 
ion intercepted X-ray luminosity is above 10 3 L© when 
q > 0.001 and above 10 4 L© when q > 0.01. These lu- 
minosities are as high as of massive main sequence stars, 
bright giants or supergiants, implying that the optical 
counterpart of a bright ULX will appear as bright as 
a massive star or giant/supergiant no matter what the 
companion really is. The observed optical spectra sug- 
gest star A has a luminosity of about 10 4 L©, and star B 
has a luminosity of about 10 3 to 10 4 £©. Thus, the 
optical emission may well be dominated by an X-ray 
"heated" companion star. However, how much of the 
companion surface is shielded by the accretion disk is 
unknown. For ULXs with Lx > 10 40 ergs s _1 , if even a 
small portion (> 10%) of the X-ray flux is not absorbed 
cnroute to the companion surface, then the companion's 
temperature and luminosity will be notably raised. 

Disk irradiation may be an important component of 
the optical emission if the X-ray luminosity is truly very 
high. Grise et al. (2007, in preparation) found significant 
variability from the optical counterpart of NGC 1313 X- 
2 on time scales from days to hours and even minutes, 
suggesting the optical emission is dominated by X-ray 
reprocessing. Both star A and star B show no detectable 
variability in the four F625W images at time scales of 
a few hours. The X-ray to optical flux ratio taking ei- 
ther star A or B as the counterpart is consistent with 
those found in neutron-star and stellar mass black hole 
X-ray binaries (the low mass X-ray binaries; LMXBs), 
where the emission is thoug ht to be dominated by repro- 
cessed X-rays on the disk (|van Paradijs fc McClintockl 
Il995t IlTu et al.ll2001[ ). However, the color of either star 
is significantly redder than the colors, (B — V)o ~ 0, 
typical of LMXBs. Assuming a standard disk geome- 
try, the disk temperature from X-ray irradiation varies 
as T(R) A oc Lx/R 2 at different radius, and the spectral 
energy distribution has a form of S v oc v~ x by integrating 
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Fig. 8. — Companion star intercepted luminosity (L c in unit of 
solar luminosity) as a function of the binary mass ratio (q, the com- 
panion to compact object mass ratio), assuming isotropic emission 
and Roche-lobe overflowing. Different lines are scaled with differ- 
ent X-ray luminosity as labeled. 

to infinite radius. Therefore, for a finite disk, the X-ray 
reprocessing at the outermost radius (i? ut) determines 
the peak of the energy distribution, which could be at- 
tributed to a simplified blackbody emission with an effec- 
tive temperature T e ff. We assume the disk size is scaled 
with the binary separation a, so the effective temperature 
of the heated disk varies as T 4 ff oc T(R out ) 4 oc Lx/a 2 . If 
the system transfers mass via Roche-lobe overflow, then 
a oc M 4 / 3 P 2 / 3 , where M is the compact object mass, P 
is the orbital period and we have assumed that the com- 
panion mass is less than that of the compact object, also 
the companion density varies as p oc P~ 2 . Thus, we find 

T ofr x Lx /4m ~ V V /6 - For most Galactic LMXBs, the 
companion is a main sequence star with a density in or- 
der of 1 gem -3 and the T e g of the heated disk is around 
25,000-30,000 K, or (B - V) ~ 0. To produce a heated 
disk with a T e g around 5,500-6,000 K as star A at a high 
X-ray luminosity level, it requires either an IMBH or an 
well evolved companion that has a low density. If the op- 
tical counterpart is star B, whose optical spectral index 
is much steeper than —1, it is unlikely that the optical 
emission is dominated by disk irradiation. 

Parenthetically, we note that the relation found by 
Ivan Paradiis fc McClintockl (|1994l ) between the visual lu- 
minosity Ly, X-ray luminosity, and orbital period of 
LMXBs has been used to make inferences about the or- 
bital periods of ULXs. This relation was based mainly on 
sources with relatively dense companions, and the corre- 
lation is inferred based on the assumption that Ly oc T 2 , 
which is no longer correct when T < 10, 000 K. There- 
fore, one should be cautious in applying this correlation 
in systems where the observed T e ff is lower than 25,000- 
30,000 K, or {B - V) is higher than 0. 

The intrinsic disk emission in the optical band should 
present in a form of S v oc u 1 ^ , which is totally different 
with the observed optical spectra. Therefore, we suggest 
that the optical emission is not dominated by the viscous 
disk. 

In sum, for a ULX system with Lx > 10 40 ergs s -1 , it 
appears unlikely that the intrinsic emission of the com- 
panion star dominates the optical emission. The com- 
panion star is likely significantly heated by X-ray illumi- 
nation. Disk irradiation could also be important. The 
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color of the optical counterpart is redder than typically 
found in systems, such as LMXBs, where a significant 
fraction of the optical light arises from an X-ray heated 
disk. The red color could be a consequence of a massive 
black hole or a low density companion. 



We thank the anonymous referee for insightful com- 
ments, which have improved the paper especially in the 
discussion section. HF thanks Manfred Pakull for helpful 
discussions. 
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